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Abstract. This article presents a phenomcnological dynainic phase transi- 
tion theory - modehng and analysis - for superfluids. As we know, although 
the time-dependent Ginzburg-Landau model has been successfully used in su- 
perconductivity, and the classical Ginzburg-Landau free energy is still poorly 
applicable to liquid helium in a quantitative sense. The study in this article is 
based on 1) a new dynamic classification scheme of phase transitions, 2) new 
timc-dcpcndcnt Ginzburg-Landau models for general equilibrium transitions, 
and 3) the general dynamic transition theory. The results in this article predict 
the existence of a unstable region H, where both solid and liquid He II states 
appear randomly depending on fluctuations and the existence of a switch point 
M on the A curve, where the transitions changes types. 



1. Introduction 

Superfluidity is a phase of matter in which " unusual" effects are observed when 
hquids, typically of helium-4 or helium-3, overcome friction by surface interaction 
when at a stage, known as the "lambda point" for helium-4. at which the liquid's 
viscosity becomes zero. Also known as a major facet in the study of quantum 
hydrodynamics, it was discovered by Pyotr Leonidovich Kapitsa, John F. Allen, 
and Don Misener in 1937 and has been described through phenomenological and 
microscopic theories. 

Atoms helium have two stable isotopes ''He and ^He. ^He consists of two elec- 
trons, two protons and two neutrons, which are six fermions. Therefore, ''He has 
an integral spin and obey the Bose-Einstein statistics. Liquid ^He, called the Bose 
liquid, displays a direct transition from the normal hquid state (liquid He I) to the 
superfluid state (hquid He H) at temperature T = 2.19K, which can be considered 
as the condensation of particles at simple quantum state. 

The main objectives of this article are 1) to establish a time-dependent Ginzburg- 
Landau model for liquid ^He, and 2) to study its dynamic phase transitions. Here- 
after, we shall present briefly the main ingredients of the study presented in this 
article. The ideas and method in this article can be used to study ^He and its 
mixture with '*Hc; this study will be reported elsewhere. 

First, in the late 1930s, Ginzburg-Landau proposed a mean held theory of con- 
tinuous phase transitions. With the successful application of the Ginzburg-Landau 



1991 Mathematics Subject Classification. 76A25, 82B, 82D, 37L. 

Key words and phrases, helium-4, dynamic phase transition, lambda point, time-dependent 
Ginzburg-Landau models, dynamic transition theory. 

The work was supported in part by the Office of Naval Research and by the National Science 
Foundation. 



2 



MA AND WANG 



theory to superconductivity, it is nature to transfer something similar to the su- 
perfluidity case, as the superfluid transitions in liquid '^He and "^He are of similar 
quantum origin as superconductivity. Unfortunately, we know that the classical 
Ginzburg-Landau free energy is poorly applicable to "^He in a quantitative sense, 
as described in by Ginzburg in [1]. 

The starting point of the modeling used in this article is to introduce a general 
principle, based on the le Chatelier principle and some general characteristics of 
pseudo-gradient flow systems. This general principle leads to a unified approach 
to derive Ginzburg-Landau type of time-dependent models for equilibrium phase 
transitions. With this general principle in our disposal, we derive some dynamic 
models for superfluid transitions. 

In addition, an important difference for the new dynamic models introduced here 
in this article from the classical ones is based on the separation of superfluid and 



normal fluid densities as described in (3.1 1 and their interactions in the Ginzburg- 
Landau free energies. 

Second, classically, phase transitions are classifled by the Ehrenfest classification 
scheme, based on the lowest derivative of the free energy that is discontinuous at the 
transition. The superfluid phase transitions has been regarded as continuous phase 
transitions, also called second-order phase transitions. However, many important 
issues are still still not clear; see among many others Ginzburg [T], Reichl [7] and 
Onuki [6^. 

One important new ingredient for the analysis is a new dynamic transition the- 
ory developed recently by the authors [3J With this theory, we derive a new 
dynamic phase transition classiflcation scheme, which classifles phase transitions 
into three categories: Type-I, Type-II and Type-III, corresponding respectively to 
the continuous, the jump and mixed transitions in the dynamic transition theory. 

The results in this article lead to two physical predictions: 1) the existence of 
the unstable region H, where both solid and liquid He H states appear randomly 
depending on fluctuations, and 2) the existence of a switch point AI, where the 
transitions, between superfluid state (liquid He H) and the normal fluid state (liq- 
uid He I), changes from first order (Type H with the dynamic classification scheme) 
to second order (Type-I). Of course, these predictions need to be verified by exper- 
iments, and it is hoped that the new model, the ideas and methods introduced in 
this article will lead to some improved understanding of superfiuidity. 

This article is organized as follows. First, the new dynamic transition theory 
is recapitulated in Section 2, and a dynamic phase transition model for ^He is 
introduced in discussed in Section 3. The dynamic transitions of the the model is 
analyzed in Section 4. Section 5 gives some physical conclusions and predictions. 

2. General Principles of Phase Transition Dynamics 

In this section, we introduce a new phase dynamic transition classification scheme 
to classify phase transitions into three categories: Type-I, Type-II and Type-III, 
corresponding mathematically continuous, jump mixed transitions, respectively. 

2.1. Dynamic transition theory. In sciences, nonlinear dissipative systems are 
generally governed by differential equations, which can be expressed in the following 
abstract form Let X and Xi be two Banach spaces, and Xi C X a compact 
and dense inclusion. In this chapter, we always consider the following nonlinear 
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evolution equations 
(2.1) 



du 



L\u + G{u,\), u{0) = if, 



where u : [0,oo) X is unknown function, and A G is the system parameter. 

Assume that Lx : Xi ^ X is a parameterized hnear completely continuous field 
depending contiguously on A G M^, which satisfies 

Lx = —A + Bx a sectorial operator, 

(2.2) A : Xi —tX a linear homeomorphism, 

Bx '■ Xi —tX a linear compact operator. 

In this case, we can define the fractional order spaces X^, for cr G M^. Then we also 
assume that G'(-, A) : X^ ^ X is C^{r > 1) bounded mapping for some < a < 1, 
depending continuously on A G M^, and 



(2.3) 



G{u,X) = oi\\u\\xJ VAg 



Hereafter we always assume the conditions (2.2) and (2.3), which represent that 
the system (2.1) has a dissipative structure. 

A state of the system (2.1 1 at A is usually referred to as a compact invariant set 
Yix- In many applications, Ea is a singular point or a periodic orbit. A state Yix of 
(2.1) is stable if Ea is an attractor; otherwise Sa is called unstable. 



Definition 2.1. We say that the system {2.1) has a phase transition from a state 
Sa at \ ~ \q j/Sa is stable on X < Xq (or on X > Xq) and is unstable on X > Xq 
(or on X < Xq). The critical parameter Xq is called a critical point. In other words, 
the phase transition corresponds to an exchange of stable states. 



Obviously, the attractor bifurcation of (2.1) is a type of transition. However, 
bifurcation and transition are two different, but related concepts. 

Let {Pj{X) G C I j G N} be the eigenvalues (counting multiplicity) of Lx, and 
assume that 



(2.4) 



(2.5) 



Re A (A) 




if A < Ao, 
if A = Ao, 
if A > Aq, 



VI < i < rn, 



yj > m + 1. 



The following theorem is a basic principle of transitions from equilibrium states, 
which provides sufficient conditions and a basic classification for transitions of non- 
linear dissipative systems. This theorem is a direct consequence of the center man- 
ifold theorems and the stable manifold theorems; we omit the proof. 



Theorem 2.1. Let the conditions (2.4-) and (2.5) hold true. Then, the system (2.1) 
must have a transition from (a, A) ~ (0, Ao), and there is a neighborhood U <Z X of 
It = such that the transition is one of the following three types: 

(1) Continuous Transition.- there exists an open and dense set Ux C U 
such that for any if dUx, the solution ux{t,Lp) of (2.1) satisfies 

lim limsup ||wA(i, <y3)||x = 0. 
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(2) Jump Transition; for any Ao < A < Ao + e with some £ > 0, there is an 
open and dense set U\ d U such that for any ip G U\, 

limsxvp \\u\{t, ip)\\x >S>0 for some 6 > is independent of X. 

t — >oo 

(3) Mixed Transition; for any Aq < A < Aq + e with some e > 0, U can be 
decomposed into two open sets Ui and U2 (U^ not necessarily connected): 
U=U^ + U^, n = 0, such that 







lim limsup iy9)||j5f 
lim sup \\u{t, </3)||x > (5 > 

t — >oo 



V<y9 e 



With this theorem in our disposal, we are in position to give a new dynamic 
classification scheme for dynamic phase transitions. 

Definition 2.1 (Dynamic Classification of Phase Transition). The phase transi- 
tions for (2.1) at X = Xq is classified using their dynamic properties: continuous, 
jump, and mixed as given in Theorem \2.1\ which are called Type-I, Type-II and 
Type-Ill respectively. 

An important aspect of the transition theory is to determine which of the three 
types of transitions given by Theorem |2.1| occurs in a specific problem. A corre- 
sponding dynamic transition has been developed recently by the authors for this 
purpose; see [3j. We refer interested readers to these references for details of the 
theory. Hereafter we recall a few related theorems in this theory used in this article. 

First, we remark here that one crucial ingredient for applications of this general 
dynamic transition theory is the reduction of (2.1 1 to the center manifold function. 
In fact, by this reduction, we know that the type of transitions for (2.11 at (0, Aq) 
is completely dictated by its reduction equation near A = Aq: 
dx 

(2.6) — = J^x + g{x, X) 



where g{x^ X) 
(2.7) g,{x,X) 



dt 

(ffi(a;,A),- 



for X e 



,5,„(a;, A)), and 



=< G(^ x,e^ + $(x, A), A), e* > 

i=l 



VI < j < m. 



Here ej and e* (1 < j < m) are the eigenvectors of L\ and L\ respectively corre- 
sponding to the eigenvalues (3j{X) as in (2.4 1, J\ is the mx m order Jordan matrix 
corresponding to the eigenvalues given by (2.4 1, and ^{x, A) is the center manifold 
function of (2.11 near Aq. 



Second, let (2.1 1 be a gradient-type equation. Under the conditions (2.4) and 
(2.5 1, in a neighborhood U C X oi u = 0, the center manifold A/^ in J7 at A = Aq 
consists of three subsets 

= ^ + ^ + 0, 

where W is the stable set, is the unstable set, and D is the hyperbolic set of 
(2.6 1. Then we have the following theorem. 



Theorem 2.2. Let (2.1) be a gradient-type equation, and the conditions (2.4-) and 
(2.5) hold true. If u = is an isolated singular point of (2.1) at X = Xq, then we 
have the following assertions: 
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(1) The transition of {2.1) at (m, A) = (0, Ao) is continuous if and only if u ^ 
is locally asymptotically stable at X = Xq, i.e., the center manifold is stable: 
iW^ = W . Moreover, (2.1) bifurcates from (0, Aq) to minimal attractors 



(2) 



consisting of singular points of (2.1) 



If the stable set of [2.1) has no interior points in M'^ 
+ D, then the transition is jump. 



I.e. 



Third, we also a dynamic transition theorem of (2.11 from a simple critical 



eigenvalue, which has been used in analyzing PVT systems and the ferromagnetic 
systems [HIS]. We refer the interested readers to these references for details for this 
theorem. 

2.2. New Ginzburg-Landau models for equilibrium phase transitions. In 

this section, we recall a new time-dependent Ginzburg-Landau theory for modeling 
equilibrium phase transitions in statistical physics. 

Consider a thermal system with a control parameter A. The classical le Chatelier 
principle amounts to saying that for a stable equilibrium state of a system S, when 
the system deviates from S by a small perturbation or fluctuation, there will be a 
resuming force to restore this system to return to the stable state S. 

By the mathematical characterization of gradient systems and the le Chatelier 
principle, for a system with thermodynamic potential Ti.{u, A), the governing equa- 
tions are essentially determined by the functional H{u, A). When the order param- 
eters • • • ,Um) are nonconserved variables, i.e., the integers 

Ui{x,t)dx — aiit) ^ constant. 



then the time-dependent equations are given by 



(2.8) 



duj 

~dt 

du , 
dn 



an 



-P,—n{u,X) + ^,{u,Vu,X) 
= (or u\dn = 0), 



for 1 < z < 



where S/Sui are the variational derivative. Pi > and 4>i satisfy 
(2.9) - - ' 



[ y"'i>i^n{u,X)dx = 0. 



The condition (2.9) is required by the Le Chatelier principle. In the concrete prob- 



lem, the terms <i>i can be determined by physical laws and (2.9 1. We remark here 



that following the le Chatelier principle, one should have an inequality constraint. 
However physical systems often obey most simplified rules, as many existing models 
for specific problems are consistent with the equality constraint here. This remark 



applies to the constraint (2.151 below as well 



When the order parameters are the number density and the system has no ma- 
terial exchange with the external, then Uj (1 < j < m) are conserved, i.e.. 



(2.10) 



Uj{x, t)dx — constant. 



This conservation law requires a continuity equation 



(2.11) 



dt 



= -V- Jj(u,A), 
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where Jj{u, A) is the flux of component Uj. In addition, Jj satisfy 

(2.12) J^ = -k,v{^i,-J2^^,), 

where /i; is the chemical potential of component ui, 

S 



(2.13) 



7Y(u, A) — (t>j{u, Vm, A), 



and 4i j{u, A ) is a function depending on the other components Ui (i j). Thus, 

dynamical equations a 

H{u, A) — (u, Vm, A) 



from (2.11 l-(2.13l we obtain the dynamical equations as follows 

" 6 



(2.14) 



du 
~di 
du 
dn 



for I < j < m, 



dAi 
dn 



\dn = 



where (3j > are constants, and 0, satisfy 



(2.15) 



n{u,X)dx ^ 0. 



When m = 1, i.e., the system is a binary system, consisting of two components 



A and B, then the term 1 



0. The above model covers the classical Cahn-Hilliard 



model. It is worth mentioning that for multi-component systems, these phij play 
an important rule in deriving good time-dependent models. 

If the order parameters (mi,-- - ,Uk) are coupled to the conserved variables 
(wfe+i, • • • , Ujn), then the dynamical equations are 



^ = ^n{u, A) + $,(m, Vu, A) 
at oui 



(2.16) 



dt 
duj 
dn 

dui 



f3,A 



-Ti{u,\) - (f>.j{u,Vu,X) 



\dn 



(or Ui\gn = 0) 







dn 



\dn 







for 1 <i<k, 

for k + I < j < m, 

for 1 < i < fc, 
for k + I < j < m. 



dn 

Here <f>i and (pj satisfy (2.9) and (2.151, respectively. 

The model (2.16) we derive here gives a general form of the governing equations 
to thermodynamic phase transitions, and will play crucial role in studying the 
dynamics of equilibrium phase transition in statistic physics. 



3. Dynamic Model for Liquid "'He 

^He was first liquidized at T = 4.215K and p ~ 1 x 10^ (Pa) by Karmerlingh 
Onnes in 1908. In 1938, P.L.Kapitza found that when the temperature T decreases 
below Tc = 2.17^, the liquid ^He will transit from normal liquid state to superfluid 
state, in which the fluid possesses zero viscosity i.e., the viscous coeflicient 77 = 0. 
The liquids with 77 = are called the superfluids, and the flow without drag is 
called the superfluidity. The superfluid transition is called A-phase transition, and 
its phase diagram is illustrated by Figure |3.1| 
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Figure 3.1. Classical PT 
phase diagram of "^He. 



3.1. Ginzburg- Landau free energy. The order parameter describing superflu- 
idity is characterized by a non vanishing complex valued function : — > C as in 
the superconductivity, originating from quantum Bose-Einstein condensation. In 
the two-fluid hydrodynamic theory of super-fluidity, the density p of ^He is given 
by 



(3.1) 



P= Ps 



where is the superfluid density and p„ the normal fluid density. The square 
is proportional to ps , and without loss of generality, we take '0 as 



(3.2) 



Based on the quantum mechanics, —ih\Iijj represents the momentum associated 
with the Bose-Einstein condensation. Hence, the free energy density contains the 
term 

where h is the Planck constant, and m the mass of atom ''He . Meanwhile the 
superfluid state does not obey the classical thermodynamic laws, which the normal 
liquid state obeys. Therefore, in the free energy density, ip satisfies the Ginzburg- 
Landau expansion 



71 



72 



2m ' ' 2 ' " ' ■ 4 
and pn has the expansion as in the free energy for PVT systems |4j . For simplicity 
we ignore the entropy, and consider the coupling action of ijj and p„, i.e., add the 
term in the free energy density. Thus the Ginzburg-Landau free energy 

for liquid ^He near the superfluid transition is given by 



(3.3) G(V,P„) = 



2m 



71 



72 



73 



Pi 2 



^.3 

3 ^" 



Ps. 4 
-jPn 



^ PnW 
P{Pn 4 



Po 



pI) 



dx. 



3.2. Dynamic model governing the superfluidity. By (2.8 1, we derive from 



(3.3 1 the following time-dependent Ginzburg-Landau equations governing the su- 
perfluidity of liquid ^He : 



(3.4) 



9V 

dpn 
dt 



kih^ 



Aljj - "flip - 72|V'PV' - IsPn-^, 



= fca Ap„ - (^1 - ppo)Pn - P2pI 



PsPn 



73 
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It is known that the following problem has a solution e H^{il) n -ffo(f^) 
standing for the density of liquid He I for any p e L^{^1): 

(3-5) dpi 



dn 



0. 



To derive the nondimensional form of (3.4 1, let 
{x,t) = ilx',Tt'), 



ml 



2po 

A2 = r(3(p°)V3 + 2p°M2 + A^oP - Ml), 



(?/;,p„) = (■(/;oV'',PoP^, + 

k2T 



where p„ is the solution of (3.5 1. 



(3.6) 



Thus, suppressing the primes, the equations (3.4) are rewritten as 
dip 



Alp + Xiip - aipnip - a2\ip\'^ip, 



dpn 

dt 



= P-Apn + AzPn - follV'P - fc2Pri - hpl 



(3.7) 



The boundary conditions associated with (3.6 1 are 
dip 



dn 



0, 



dpn 

dn 



= 



on 9ri. 



When the pressure p is independent of x £ fi, then the problem (3.6) and (3.7) 
can be approximatively replaced by the following systems of ordinary differential 
equations for superfluid transitions: 
(hp_ 
Hi 

dpn 

dt 



(3.8) 



AiV' - aipnip - a2\ip\^ip, 
= A2P„ - hi\ip\'^ - b2pi - hpl- 



By multiplying ip* to both sides of the first equation of (3.8) and by (3.2 1, the 

o: 

= AiPs 



aiPnPs - a2Ps, 
= A2P„ - bips - b2pl ~ b3p 



equations (3.8) are reduced to 

dps 
dt 

(3.9) dp, 

dt 

{Ps{0),Pn{0)) = {xo,yo)- 

for Ps >0 and xq > 0. 

We need to explain the physical properties of the coefficients in ( 3.4 1 and ( 3.6 ) . It 
is known that the coefficients 7i (1 < « < 3) and pj (0 < j < 3) depend continuously 
on the temperature T and the pressure p: 

lz^li{T,p), p-j = Pj{T,p) VI < z < 3, 0<j<3. 
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From the both mathematical and physical points of view, the following conditions 
are required: 

(3.10) 72 >0, 73 >0, fi3>0 \fT,p. 

In addition, by the Landau mean field theory we have 

r >0 if{T,p)eAu 



(3.11) 



(3.12) 



7i 



Ml 



<0 if(T,p)eA2, 

>0 if(T,p)eBi, 
<0 if(T,p)eB2, 



where Ai, Bi (i = 1, 2) are connected open sets such that Ai -\- A2 — Bi + -B2 — I^+i 
and Ai n ^2, ^1 n B2 are two simple curves in M^; see Figure 4.1 (a) and (b). In 
particular, in the PT-plane, 



(3.13) 



97 



dT 



i>0, 



dp 



i>0, ?^^>0 



9/i 



dT 



dp 



i>0. 



By (3.101 the following nondimensional parameters are positive 
(3.14) 



ai > 0, 02 > 0, 61 > 0, 63 > 0. 



By (3.11l-(3.13l, the following two critical parameter equations 

Ai = Ai(T,p) = 0, 
A2 = A2(T,p) = 0, 



(3.15) 



4.1 



give two simple curves li and I2 respectively in the PT-plane M^; see Figure 
(c) and (d). 

The parameters /i2 and 62 depend on the physical properties of the atom He, 
and satisfy the following relations: 

62(7",^) = t/9o(M2 + 3p°/i3) < iff psoi > /O; at A2(T,p) = 0, 

(3.16) 

b2{T,p) = r/9o(/^2 + 3p„/i3) > iff psoi < Pi at \2{T,p) = 0, 



where psoi and pi are the densities of solid and liquid, and the solution of (3.5) 
representing the liquid density. These relations in (3.16) can be deduced by the 
dynamic transition theorem of (2.1 ) from a simple critical eigenvalue, Theorem A. 2 
in Ma and Wang [5]. 

4. Dynamic phase transition for liquid "'He 

In order to illustrate the main ideas, we discuss only the case where the pressure 
p is independent of a; € fi, i.e., we only consider the equations (3.8 1 and (3.9); the 
general case can be studied in the same fashion and will be reported elsewhere. 



4.1. PT-phase diagram. Based on physical experiments together with (3.11)- 
(3.13), the curves of ^i{T,p) = and pi{T,p) = in the PT-plane are given by 
Figure 4.1 'a) and (b) respectively. By the formulas 

Ai(T,p) = -t(7i(T,p) + 73P^), 

^2{T,p) = -t{pi{T,p) - pop-2plp2 - 3p"V3), 

together with (3.10) and ( 3.16[ ), the curves li and I2 given by (3.15) in the PT-plane 
are illustrated in Figure |4.1[c) and (d) . 
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Figure 4.1. (a) The curve of 71 = 0, (b) the curve of /ii = 0, (c) 
the curve of Ai = 0, (d) the curve of A2 = 0. 



Figure 4.2. ACF is the curve of 
A2 = 0, BCD is the curve of Ai = 
0. 



Di = {iT,p)eRl I Ai(r,p) >0}, 
D2^{{T,p)eRl I Ai(r,p) <0}, 
E,^{iT,p)eRl I A2(r,p)>0}, 
E2^{iT,p)eRi I A2(r,p) <0}. 



Let the curve of Xi{T,p) = intersect with the curve of X2{T,p) = at a point C; 
see Figure |4.2[ 
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When Ai crosses the curve segment CB to enter into Di from D2 (see Figure 



4.2 1, we have A2 < and Ai satisfies 



Ai(T,p) 




if (r,p) eD2, 
if (r,p) e CB, 
if (r,p) ei^i. 



In this case, by Theorem 2.1 the equations ( |3.9| have a phase change which de- 
scribes the transition from hquid He-I to hquid He-II, and the second equation of 
(3.9 1 can be equivalently rewritten as 

dt 

where pn — Pn — Pn^ Pn satisfies the equation 

Mpn - hpl - bspl = bips 



(4.1) 



= (A2 - 2b2pl - Sh{plf)K ~ {b2 + ?>hpl)pi - hpl 



and ps > is the transition solution of (3.9). Since p* is the density deviation of 
normal liquid, by (3.1 1 we have 

Pn^ P~ Pn- Ps- 

On the other hand, (~ p) represents the density of liq uid He I. Thus, p* ~ 
—Ps < 0. From the PT-phase diagram of ^He (Figure 3.11, we see that pg^i > Pi 
for "^He near T = 0, therefore by (3.16), 62 < 0. Hence we derive that 

(4.2) A2-2fe2P;-363(p;)' <0, 



for {T,p) in the region Di \ H, as shown in Figure 4.2 



It follows from (4.1 1 and (4.2) that when {T,p) is in the region Di \H, the liquid 
^He is in the superfluid state. 



When A2 crosses the curve segment CA, as shown in Figure |4.2| to enter into 
El from E2, then Ai < and 

<0 if iT,p) e E2, 
= if (T,p) e CA, 
> if {T,p) e El. 



In this case, the equations (3.9) characterize the liquid-solid phase transition, and 
the first equation of (3.9 1 is equivalently expressed as 

dps 
dt 

where p„ > is the transition state of solid ''He . By (3.14) we have 
(4.4) 



(4.3) 



= (Ai - aipn)ps 



a2Ps 



Xi(T,p) - aipn < 0, 



for any (T,p) in the region Ei \ H , as shown in Figure 4.2 



From (4.3) and (4.4) we can derive the conclusion that as {T,p) in Ei\H,ps = 
is stable, i.e., *He is in the solid state. 

However, the shadowed region H in Figure [4!2] is an unstable domain for the solid 
and liquid He II states, where any of these two phases may appear depending on the 
random fiuctuations. Thus, from the discussion above, we can derive the theoretical 
PT-phase diagram given by Figure [473| based on equations (3.9). In comparison 
with the experimental PT-phase diagram (Figure 3.1 1, a sight difference in Figure 
|4.3|is that there exists an unstable region H, where the solid phase and the He II 
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Figure 4.3. Theoretical PT- 
phase diagram: H is an un- 
stable region where both solid 
and liquid He II states appear 
randomly depending on fluctu- 
ations. The point M is a point 
where the transitions, between 
superfluid state (liquid He II) 
and the normal fluid state (liq- 
uid He I), changes from Type 
II to Type-I; namely, the tran- 
sition crossing CM is Type- 
II and the transition crossing 
MB is Type-I. 
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phase are possible to occur. This unstable region in Figure 4.3 corresponding to 
the coexistence curve CE in Figure 3.1 is a theoretical prediction, which need to 
be verified by experiments. 

4.2. States in unstable region. We consider the dynamical properties of transi- 
tions for (3.9) in the unstable region. It is clear that at point C = (TcPc), 

(4.5) Xi{Tc,pc)^0, HTc,Pc)^0, 

and the unstable region H satisfies that 

H =. {{T,p) e Rl I Ai(r,p) > 0, A2(T,p) > 0}. 



To study the structure of flows of (3.9) for {T,p) e it is necessary to consider 
the equations (3.9) at the point C = (TcPc), and by (4.5 1 which are given by 

dps 
dt 
dpn 
dt 



(4.6) 



-aiPnPs - a2Ps 



Due to (3.16) and p^oi > pi, we have 
(4.7) 



62 < for {T,p) C H. 



Under the condition (4.7), equations (4.6) have the following two steady state so- 
lutions: 

Zi = {ps,pn) = (0,1621/63), 

Z2 = {Ps,Pn) = ( —a, -a] , 



\b2\ 

263 



4ai6i63/a2|62p - !)■ 

By direct computation, we can prove that the eigenvalues of the Jacobian matrices 
of (4.6) at Zi and Z2 are negative. Hence, Zi and Z2 are stable equilibrium points 
of (4.6). Physically, Zi stands for solid state, and Z2 for superfluid state. The 
topological structure of (4.6) is schematically illustrated by Figure 4.4 ^a), the two 
regions Ri and R2 divided by curve AO in Figure [4!4] (b) are the basins of attraction 
of Zi and Z2 respectively. 
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We note that in H, Ai and A2 is small, i.e., 
0<Ai(T,p), A2(T,p)«l, 



for (T,p) e H, 



and (3.9 1 can be consider as a perturbed system of ( |4.6| . 

Thus, for (T, p) G H the system (3.9 1 have four steady state solutions Zi — 
Z{T,p) (1 < i < 4) such that 

lim (Zi(T,p),Z2(T,p),Z3(r,p),^4(r,p)) = (^1,^2,0,0), 

and Zi and Z2 are stable, representing solid state and liquid He-I I sta te respectively, 
and Zi are two saddle points. Th e topological structure of (3.9 1 for {T,p) e H 
is schematically shown in Figure 4.4 c) , an d the basins of attraction of Zi and Z2 
are i?i and R2 as illustrated by Figure 4.4 d). 





Figure 4.4. 



In summary, with the above analysis and the dynamic transition theory, we 
arrive at the following theorem: 

Theorem 4.1. There exist Jour regions Ei,Di and H in the PT-plane (see Figure 
4- -3), which are defined by 



E, = {iT,p)e 
Di^{{T,p)e 
^2 = {(T,p)e 
H = {{T,p)e. 

such that the following conclusions hold true: 

(1) If{T,p) G El, the phase of^He is in solid state. 

(2) If(T,p) e Di, the phase is in superfiuid state. 

(3) If(T,p) G D2, the phase is in the normal fluid state 



ll I Ai(T,p) <0, A2(r,p) >0}, 
ll I Ai(r,p) >0, A2(T,p) <0}, 
ll I Ai(r,p) <0, A2(T,p) <0}, 
I I Ai(r,p) >0, A2(T,p) >0}, 
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(4) If {T,p) G H, there are two regions Ri and R2 in the state space {ps,Pn) 
suc h th at, under a fluctuation which is described by the initial value (xq^i/q) 



3.9), if (xq, uq) G Ri then the phase is in solid state, and if (xq, i/q) £ R2 



then it is in superfluid state. 

In general, the observed superfluid transitions of ''He are of the second order, 
i.e., Type-I (continuous) with the dynamic classification. But, from (3.9 1 we can 
prove the following theorem, which shows that for a higher pressure the superfluid 
transitions may be of zeroth order, i.e., the Type-II (jump) transition with the 
dynamic classification scheme. 



Theorem 4.2. Let {Tq,pq) satisfy that \i{To,Pq) = 0, A2(To,po) < 0. Then, pjp 
have a superfluid transition at (To,po) from D2 to Di. In particular, the following 
assertions hold true: 

(1) Let 

^^TrT~"2 at {T,p) ^ (To,pq). 
IA2I 

Then if A < the superfluid transition is Type-I, and if A > the superfluid 
transition is Type-II. 

(2) The equation 



(4.8) 



determines a point M on the CB coexistence curve in Figure \4_ 
the transition changes type from Type-II to Type-I. 

Proof. Step 1. By the assumption, it is clear that 

r <0 if {T,p)£D2, 
Ai(T,p)= J =0 if (T,p) = (To,po), 
I >0 if {T,p)£Dr, 



whe 



where Di,D2 are as in Figure 4.1 'c). Hence, by Theorem 2.1 (3.9) have a tran- 



sition at {Tq,pq). By X2{To,Po) < 0, the first eigenvalue of (3.9 1 is simple, and its 
eigenvector is given by e — (e^, e„) — (1,0). 

The reduced equation of (3.9 1 on the center manifold reads 

dx 



dt 



— Xix — aixh{x) — a2X , x > 0, 



where h(x) is the center manifold function. By (3.9), h can be expressed as 



h{x) 



A2 



X + o{x'^), X > 0. 



Thus the reduced equation of (3.9) is given by 
dx 



(4.9) 



Xix + Ax^ + 0(2:^) 



for X > 0, 



where A is as in (4.8). The theorem follows from (4.9 1 and Theorem 2.2 



Step 2. By the nondimensional form, we see that 
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where 73 is the coupled coefRcient of p„ and pa- Physically, 73 is small in comparison 
with 72; namely 

< aihi < 02. 

On the other hand, we know that 

A2(To,po)^0 as (To,Po) ^ (7c,Pc), 



where C — (TcPc) is as in (4.51. Therefore we deduce that there exists a pressure 

p*{b < p* < pc) such that 

Qibi a [ if 0<po<P*, 

|A2(7o,po)| ^[ >0 if p*<po<pc• 
Thus, when the transition pressure po is below some value po < P* the superfluid 
transition is of the second order, i.e., is continuous with the dynamic classification 
scheme, and when po < p* the superfluid transition is the first order, i.e., is jump 
in the dynamic classification scheme. The proof is complete. □ 

5. Physical Conclusions and Predictions 

As we know, the classical phase transition is illustrated by Figure |3.1[ In this 
phase diagram, the coexistence curve CE separates the solid state and the super- 
fluid state (liquid He II), and the curve CB is the coexistence curve between the 
superfluid state (liquid He II) and the normal fluid state (liquid He I), and the 
critical point is the triple-point. It is considered by the classical theory that the 
transition crossing CB is second order in the Ehrenfest sense (Type-I with the 
dynamic classification scheme). 



However, in the phase transition diagram Figure |4.3| derived based on Theo- 
rems 



4.1 and |4.2| there is a unstable region H, where both solid and liquid He II 
states appear randomly depending on fluctuations. The point M is a switch point 
where the transitions, between superfluid state (liquid He II) and the normal fluid 
state (liquid He I), changes from first order (Type II with the dynamic classifi- 
cation scheme) to second order (Type-I); namely, the transition crossing CM is 
second order (Type-II) and the transition crossing M B is first order (Type-I) . 

In summary, the results in this article predict the existence of the unstable region 
H and the existence of the switch point M. It is hoped these predictions can be 
verified by experiments. 
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